Classical solutions from quantum regime for barotropic FRW model 
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The quantization of gravity coupled to barotropic perfect fluid as ma- 
tter field and cosmological constant is made and the wave function can be 
determined for any k in the FRW minisuperspace model. The meaning of 
the existence of the classical solution is discussed in the WKB semiclassical 
approximation . 
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I. INTRODUCTION 

It is a general belief that the existence of an initial singularity can be remove through 
the employment of a quantum theory of gravity. However, there is no consistent theory of 
gravity until now, and in this sense the problem of the initial singularity remains of actuality. 

It is well known that is possible to construct a quantum model for the universe as a 
whole, through the Wheeler-DeWitt (WDW) equation, based in the ADM decomposition 
of the gravity sector, which leads to a hamiltonian approach of general relativity, from 
which a canonical quantization procedure can be applied (Halliwel, 1991). Moreover, in the 
hamiltonian formalism, the notion of time is lost (Isham, 1992), but there are some proposals 
by which this notion of time can be recovered (Schutz, 1970, 1971) coupling of the gravity 
sector to a perfect fluid. In this scheme, called Schutz's formalism, a quantization procedure 
is possible, and the canonical momentum associated with the perfect fluid appear linearly in 
the WDW equation, permiting to rewrite this equation in the form a Schrodinger equation 
with a time coordinate associated with the matter field. 

Recently, a convincing quantum gravity origin for inflation should however be proved. 
According to Bojowald (Bojowald, 2002) a non-perturbative approach would produce the 
most reliable answer as to whether or not inflation can be derived from quantum gravity. 
At the less ambitious level of minisuperspace quantum cosmology the generality of inflation 
as well as the quantum creation of the universe should be treated by analyzing the Wheeler- 
DeWitt equation and its physical solutions, the so-called wave functions of the universe 
(Gibbons and Grishchuk, 1989). 

In this way, the behaviour of the scale factor may be determined in two different forms 
from the quantum regime: when the wave function is time dependent, we can calculate 
the expectation value of the scale factor, in the spirit of the many worlds interpretation of 
quantum mechanics ( Tipler, 1986) 

J "W(A)**(A,t)A(t)tt(A,t)dA 

/ ~W(A)tf*(A,t)tf(A,t)dA ' 1 ' 
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where W(A) is a weight function that normalizes the expectation value. The other way 
is to apply the WKB semiclassical approximation, in which case the system follows a real 
trayectory given by the equation 

where the index q designates one of the degree of freedom of the system, and <3> is the phase 
of the wave function when is written as 

$ = We i$ , (3) 

the functions W and <3> are real functions. We follow this last procedure to obtain the 
classical behaviour for the scale factor, using the FRW model with a barotropic perfect fluid 
and cosmo logical constant. 

The work is organized as follow. In next section we describe the quantum model with 
the solution for any n case in the minisuperspace, considering a barotropic perfect fluid as 
matter field including the cosmological constant. In section III, we present the classical 
evolution for the scale factor derived following the semiclassical WKB procedure. Section 
IV is devoted to conclusions. 



II. QUANTUM MODEL 

The total lagrangian for a barotropic perfect fluid coupled to gravity using the FRW 
geometry with the classical cosmological constant term is given by 

^total £geom ^matter- (4) 

namely 

2 



I 777 r nA1 6A 2 d 2 A 6A /dA\ 6A 2 dAdN . . . ^ 

"""" = ^\K-2K] = -—^--[-) +w --- 6*NA - 2NAA 

d/-6A^ + 6A/dAY_ feNA _ 2NAA3i (5) 



dt V N J N V dt J 
and the lagrangian matter density (Ryan, 1972; Pazos, 2000) 

Cmatter = V^^g [l67rGNp { ( 7 + 1) (l + g km U k U m ) k - 7 (l + g km U k U m ) ~* } 

-16ttGp( 7 + l)U m N m ] . (6) 
These results are obtained by employing the ADM form of the FRW metric 

dr 2 



ds 2 = — N dt" + A 



1 — KT' 



+ r 2 (d^ 2 + sm 2 9d<t> 2 ) 



(7) 



where N is the lapse function, A is the scale factor of the model, and k is the curvature 
index of the universe ( k — 0, +1, —1 plane, close and open, respectively). 
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We also make use of the usual perfect fluid energy-momentum tensor 



V = pg M , + (p + p) U M U, , (8) 

where p, p, are the pressure, energy density and the four-velocity of the cosmological 
fluid, respectively; and using the barotropic relationship p = 7p, 7 = constant, we have the 
solution for the energy density as a function of the scale factor of the FRW universe, in the 
usual way 

P ~ A3(7+l) ' ^ 

where M 7 is an integration constant. In C ma tter we also choose the comoving fluid (three- 
velocity Uk = 0), and the gauge N k = 0, obtaining 

Cmatter = 16ttGNM 7 A^. (10) 

Thus, the total lagrangian density has the following form 

Ctot = ± (^) + ^ (f ) 2 - 6-NA - 2NAA 3 + 167rGNM 7 A~ 37 . (11) 

Following the well-known procedure to get the canonical Hamiltonian function, we define 
the canonical momentum conjugate to the generalized coordinate A (scale factor) as 

n A .|, ( 12 > 



L = n A A - NTi — n A A - N 



TT 2 

— ^ + 6kA + 2AA 3 - 16ttGM 7 A- 37 
24A 7 



(13) 



where 



n 



TT 2 

24A 



+ 6kA + 2AA 3 - 16ttGM^A" 37 . 



(14) 



Performing the variation of (13) with respect to N, ^ = 0, implies the well-known result 
Ti, — 0. Imposing the quantization condition and applying this hamiltonian to the wave 
function we obtain the WDW equation in the minisuperspace 



H\V) 



1 



24A 



dA 2 



+ 144kA 2 + 48AA 4 - 384ttGM^A- 37+1 



= 0. 



(15) 



Notice that in principle the order ambiguity in equation (15) should be taken into account. 
This is quite a difficult problem to be treated in all its generality, since the hamiltonian 
operator in (15) must be written in a very general form in order to take into account all 
possible order, but, at least in the minimal case in which (Hartle and Hawking, 1983) (there 
are other possibilities depending on different considerations on the operators, see the works 
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of Christodoulakis and Zanelli, (Christodoulakis and Zanelli, 1984) or Lidsey and Moniz 
(Lidsey and Moniz, 2000) 

. ,d 2 ^ . , , n d . n d^ . i /d 2 ^ . ,d#\ 
A " dA* ~ dA dA = (dJF - PA dA J ' (16) 

where the real parameter p measures the ambiguity in the factor ordering. Therefore, the 
Wheeler-DeWitt equation can be written as follows 

d 2 ^ dty 

- A lA^"lA- V ^ = °- < 17 > 

where 

V(A) = -48AA 5 + 3847rGM 7 A~ 37+2 - 144reA 3 . (18) 

In the following we discuss some of the quantum solutions of equation (17) for particular 
values of the 7 parameter and the parameter p of factor ordering. 

A. Inflationary scenario 

In the inflationary era, including the cosmological constant A, we choose 7 = — 1. Notice 
that the WDW equation (17) can be written as 

- P^r + 144A 3 f m 2 A 2 - kW = 0, (19) 
dA 2 dA v > 

where 

m 2 = ~ + ^ttCM.! . (20) 
The latter relationship leads to three possible cases 

1. m 2 > 1 

The differential equation for this subcase is 

Atf" ! - Q*'_! + 144kA 3 (m 2 A 2 - nj *_i = (21) 

which, after the substitutions v = m 2 A 2 — k, = v 1 / 2 y(v), and z = ^v 3 / 2 , turns 
for p = 1 into a Bessel equation with the general solution 

4 r o . o 1 3/2 



*_i = fm 2 A 2 - k) 2 [a Ji (z) + b J_i (z)l , where z = — fm 2 A 2 - k] (22) 

V / L 3 3 J xx^2 L J 

where ao and bo are superposition constants. 



1 In general m 2 is not a positive constant, see the corresponding definition (20). 
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2. m 2 < 0, 

The differential equation for this subcase is 

-Atf'^ + p*'_! + 144/cA 3 (|m 2 |A 2 + k) = 0. 

For p = 1 we have a modified Bessel equation with the general solution 

4 



= (|m 2 |A 2 + K) 2 [aili(z) +biKi(z) 



where 



z = 



nr 



m 2 |A 2 + k 



(23) 



3/2 



(24) 



where ai and bi are superposition constants. 
3. For m 2 = 0, the differential equation for this situation is 

-A^f'U + p*'_! + 144kA 3 ^_i = 

which for k — 1, has as solution the modified Bessel functions of order v 

= A 2 ^ fA I,(6A 2 ) + B K,(6A 2 ) 



i+p 
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(25) 



(26) 



where A and B are superposition constants, while for k — — 1, the solution become 
to be the ordinary Bessel functions 



*_i = A 2u [AiJ„(6A ) + BiYj,(6A ) 
where Ai and Bi are superposition constants. 

B. Inflationary like scenario, 7 = — 

In this particular case, the WDW equation becomes 
d 2 ^ d^ 



(27) 



A 



- p— - 48A 3 (AA 2 - g) * = 0, with g = 87rGM_ 1/3 - 3k 



dA 2 *dA 



(28) 



which, after the substitutions v = AA 2 — g, ^ = v 1 / 2 y(v), and z = ^^v 3 / 2 , yields for p = 1 
a Bessel equation with the general solution 



*(A) = AA 2 - g 



2 ' ' k' I1 Ul^rAA a - g rUc 1 K 1 



1/2 



3 y 3A 

where Co and Ci are superposition constants. 



lv'3 t v2 



3 V 3A 



AA^-g 



3/2 



• (29) 
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C. Dust case (7 = 0, k = 0) 
For this case the WDW equation is 

A S ~ p dT + 48A2 (~ aa3 + 87rGM °) * = °' ( 3 °) 

Making the transformations z = ^y/3AA 3 and ^ = e~5 z u(z), one gets for u(z) the confluent 
hypergeometric equation 

d 2 u , . du n . . 

z d^- z >s- mu = (31) 

where m = ^=E _ 16 ^| l0 and 7 = 2 -y^. The linear independent solutions are (Gradshteyn 
and Ryzhik, 1980) 

Ui(z) = 1 F 1 (m, 7 ;z) (32) 
u 2 (z) =z 1 -^ 1 F 1 (m- 7 + 1,2-7; z) (33) 

where 1 F 1 is the degenerate hypergeometric function. 
Thus, the exact solution for \& becomes 

tt(A) = e^ z [Aoux(z) + B u 2 (z)], (34) 

where A and B are superposition constants, ui and u 2 are the functions given in (32) and 
(33), respectively. 



D. Stiff" fluid case (7 = 1, k = 0) 

The WDW equation 

A 3T? - PTT + 48 (-AA 5 + SttGMxA- 1 ) = 0, (35) 
dA 2 dA v > 

has the following exact solution 

$(A)=A^Z„^A 3 J, with ^ = ^(^) 2 - 3847TGML (36) 

The exact expression of the wave function depends on the sign of the cosmological constant 
and the order v: 

• A > and v real, the functions Z„ become the modified Bessel functions, either l u or 
K u , depending on the boundary conditions. 

• A < and v real, the functions Z„ turn into the ordinary Bessel functions, either J„ 
or Y„, depending on the boundary conditions. 

• A > and v pure imaginary, the functions Z u become the modified Bessel functions 
of pure imaginary order (Dunster, 1990), either \ v or K u , depending on the boundary 
conditions. 

• A < and v pure imaginary, the functions Z„ turn into the ordinary Bessel functions 
of pure imaginary order (Dunster, 1990), either J u or Y u , depending on the boundary 
conditions. 
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III. THE CLASSICAL BEHAVIOUR FROM WKB REGIME 



Interesting results can be obtained at the level of WKB method if one performs the trans- 
formation n A — > ■ Then, (14) becomes the Einstein-Hamilton- Jacobi equation, where $ 
is the superpotential function that is related to the physical potential under consideration. 

Introducing this ansatz in (14) we get 



H 



24A 
thus, we obtain 



^ + 144/tA 2 + 48AA 4 - 384ttGM 7 A- 37+1 



0, (37) 



d$ / 

— = V /-48AA 4 + 384vrGM 7 A-^+ 1 - 144kA 2 (38) 



the superpotential <3> has the following form 



$ = ± / ^/-48AA 4 + 3847rGM 7 A-^ +1 - 144/tA 2 dA, (39) 

This integral can be solved for particular cases of the 7 parameter. 

Now, employing the equations (2,12) and (38) we will obtain the evolution for the scale 
factor. The classical equation of motion is 



1 9 A rl A / 1 8 

— ^ = U A = VA^-AA* + %GM,A-^ - K , (40) 

in term of the "cosmic time" r defined by dr = N(t)dt, the equation (40) is read how (when 
we choose the gauge N(t) = l, this cosmic time is the physical time t). 

dA , . 

dr = . 41 

^/-|AA 2 + §ttGM 7 A- 3 7-i - k 



A. Inflationary scenario, 7 



The equation (41) for the inflation regime is written as 
rA dx 1 



'° ^/(-|A + §7rGM 7 )x 2 - k 
where ai = — |A + |7rGM 7 and a\ = 



In 



A + v /A 2 -a 2 



(42) 



The inverse of the equation (42) us give the following structure for the scale factor A(r) 



AM = - 



eA /aT(r-To) _|_ & 2 e -y/Ei(T-To) 



(43) 



The scale factor A(r) will have a inflationary behaviour only if the parameter ai 3> 1, then 
the cosmological constant have the restriction value A < 3 (|7rGMi — l). If the parameter 
< a! < 1, the scale factor vanishes very fast. 

For other scenarios and particular values in k and A = 0, we have the following: 
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1. For the dust case and k — 0, we obtain the well-known solution A oc t 2 / 3 for N=l. 

2. For 7 = 1 and k — 0, the scale factor have the traditional solution A oc t 1 / 3 for N=l. 

3. For radiation case and k = — 1, the behaviour is 



A(r) = y4(r-r ) 2 --^GMi. (44) 

For the case of k = 0, A oc t 1 / 2 for N=l. 

On the other hand, the master equation (41) can be solved for any 7 if the cosmological 
constant vanish, defining dr = A 37+2 dT, having the following general solution in the time T 



/ 3 \ 3 7 +l 

A ' T » = (sgs;) 



{i I rGM 7 (3 7 + l)(T-T„)J +k 



3 7 +l 

(45) 



The classical solutions obtained by this method perhaps coincide with these obtained 
solving the Einstein field equation, but is not possible to know if the quantum universe 
remains quantum forever for some case, due to the fact that the WDW equation is not time 
dependent. 

IV. CONCLUSIONS 

We discussed from the point of view of simple models of minisuperspace quantum cos- 
mology. We found that the wave functions of the WDW equations in question are mostly 
Bessel functions. For particular values of the parameter p that measures the ambiguity 
of the factor ordering problem we obtained analytical results for a some scenarios of the 
universe possessing a cosmological constant. Only in the case of stiff fluid, the cosmological 
constant can be positive or negative. In the other cases, it can be only positive for physical 
reasons related to the scale factor or the potential under consideration. It is well known that 
the Wheeler-DeWitt cosmological equation is not an evolution equation and therefore the 
associated quantum states do not evolve in time. A possible way out of this difficulty could 
be to connect some parameters of the 'quantum' WDW solutions with classical Einstein 
ones by phenomenological restrictions imposed on the superpotential function as we did in 
this work. Using this method, we find the classical behaviour for the scale factor, that are 
analogous to those found solving the Einstein equations. 
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